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Background 

Myelodysplastic syndrome (MDS) is a preleukemic condition that transforms into acute 
myeloid leukemia. However, the genetic events underlying this transformation remain 
poorly understood. Aberrant DNA methylation may play a causative role in the disease 
and its prognosis. Thus, v\/e compared the DNA methylation profiles in refractory anemia 
\N\t\\ excess blast (RAEB) to those in refractory cytopenia Wxth multilineage dysplasia 
(RCMD). 

Methods 

Bone marrov\/ samples v\/ere collected from 20 patients W\t\\ primary MDS (9 W\t\\ RAEB 
and 1 1 W\t\\ RCMD), and peripheral blood samples w/ere collected from 4 healthy controls. 
These samples wieve assessed using a commercial v\/hole genome-v\/ide methylation 
assay. Methylation-specific polymerase chain reaction (PGR) was used to detect the meth- 
ylation of candidate gene promoters in RAEB and RCMD. 

Results 

Microarray data revealed significant hypermethylation in 69 genes within RAEB but not 
RCMD. Candidate genes were mapped to 5 different networks, and network 1 had the 
highest score due to its involvement in gene expression, cancer, and cell cycle. Five genes 
{GSTM5, BIK, CENPH, RERG, and ANGPTL2) were associated with malignant disease 
progression. Among them, the methylated promoter pairs of GSTM5 (55.5% and 20%), 
BIK (20% and 0%), and ANGPTL2 (44.4% and 1 0%) were observed more frequently in 
RAEB. 

Conclusion 

DNA methylation oiGSTMB, BIK, SindANGPTL2 may induce epigenetic silencing and con- 
tribute to the increasing blasts and resulting MDS progression; however, the functions 
of these genes were not determined. Further study focusing on epigenetic silencing using 
various detection modalities is required. 
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INTRODUCTION 



Myelodysplastic syndrome (MDS) is a preleukemic con- 
dition that transforms into acute myeloid leukemia (AML) 
in approximately 30-40% of cases. MDS-associated leukemia 
is rarely cured by chemotherapy, and patient prognosis is 



typically poor [1]. MDS represents an excellent model of 
leukemic development with progressive increases in blastic 
bone marrow involvement; however, the genetic events un- 
derlying this evolution remain poorly understood. 

MDS is separated into subtypes based on bone marrow 
blast and dysplasia. Among the MDS subtypes, refractory 
anemia with excess blast (RAEB) is the most common 
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pre -leukemic type and is associated with poor prognosis. 
Thus, the focus of the present study is on the bone marrow 
blast subtype of MDS, specifically on different DNA methyl- 
ation profiles between the subtypes morphologically affected. 

A whole genome-wide DNA microarray technique was 
used in the present study to assess the DNA methylation 
profiles between RAEB and refractory cytopenia with multi- 
lineage dysplasia (RCMD). Genes displaying different DNA 
methylation patterns were identified using Ingenuity 
Pathway Analysis (IPA) tools. In this study, we aimed at 
identifying genes that were hypermethylated in RAEB but 
not in RCMD. 



MATERIALS AND METHODS 



1. Sample collection 

This study included 20 patients with primary MDS. As 
per the WHO criteria, 11 patients had RCMD and 9 had 
RAEB. The study was approved by the ethics committee 
of Dong-A University Hospital, and informed consent was 
obtained from each participant. Following consent, a bone 
marrow aspirate was collected from each patient, and periph- 
eral blood samples were harvested from 4 healthy volunteers. 
Density gradient centrifugation was used to separate mono- 
nuclear cells from bone marrow or peripheral blood 
specimens. Genomic DNA was extracted using the QIAmp 
DNA MiniKit (Qiagen, Valencia, CA, USA). 

2. DNA methylation array 

Methylation levels were detected using the Infinium 
Methylation Assay (lUumina, San Diego, CA, USA) according 
to the manufacturer's protocol. Briefly, DNA was treated 
with sodium bisulfate from an lUumina-specified kit (Zymo 
Research, Orange, CA, USA), which converted unmethylated 
cytosine to uracil and left methylated cytosine unchanged. 
For each CpG site, 2 pairs of probes were designed that 
corresponded to either the methylated or unmethylated state 
of the CpG site. Polymerase chain reaction (PGR) templates 
were generated via allele -specific extension and ligation. 
These templates were PGR amplified using fluorescent-la- 
beled common primers. The resulting PGR products hybri- 



dized to a bead array at sites with complementary sequences, 
and the fluorescent label of the hybridized targets denoted 
a methylated or unmethylated state for a given locus. The 
methylation status of the interrogated CpG site-the p-val- 
ue-was represented by the ratio of fluorescent signal from 
1 allele relative to the sum of both methylated and un- 
methylated alleles. The p-value provided a continuous meas- 
ure of DNA methylation at a CpG site and ranged from 
0 (no methylation) to 1 (complete methylation). The overall 
average p-value was provided for all 27,578 CpG sites found 
in the 14,495 genes analyzed from a particular sample. In 
this study, the methylation status of CpG sites in both groups 
was assessed based on the median p-value (0.5). The p-value 
of hypermethylated genes was greater than 0.5 in the RAEB 
group samples and less than 0.5 in the RCMD group samples. 
The significance of the difference in methylation status be- 
tween RAEB and RCMD was assessed by 4 different methods. 
The Mann- Whitney and Kmskal-Wallis tests were used for 
nonparametric analysis, and the t-test and ANOVA were 
used for parametric analysis. 

3. IPA 

IPA tools were used to conduct a bioinformatics evaluation 
of methylated candidates [2]. The core bioinformatics analy- 
ses were performed using the Pathway Studio V.5.0 (Ariadne 
Genomics, Rockville, MD, USA). IPA contains a knowledge 
repository of oncological networks and biological relation- 
ships that have been encoded systematically based on more 
than 200,000 original peer-reviewed articles. Details and 
technical requirements are available from either the IPA 
Website (www.ingenuity.com) or the supplemental methods 
in the online appendix. A graphical representation shows 
the IPA-generated hypothetical networks based on the mo- 
lecular relationships, interactions, and pathway associations 
among methylated candidate genes. In this study, IPA tools 
were used to analyze significant differences in gene methyl- 
ation and select genes with different methylation levels in 
RAEB and RCMD. 

4. Methylation-specific PCR (MSP) analysis 

MSP was used to validate candidate methylated genes 
in samples obtained from the same set of 20 MDS patients 



Table 1. Methylated allele p 


rimer sequences and optimized annealing temperatures for methylation-specific PCR. 




Gene 


Sequence 


CpG sites 


Length (bp) 


Annealing temp CC) 


BIK 


F: 5' TTTGTTCGTTTCGTTAGGC 3' 


3 


541 


58 




R: 5' TTTATCCCGTACGTCGCTAT 3' 








ANGPTL2 


F: 5' ATGAGGTTTTTACGAAGGGC 3' 


2 


368 


60 




R: 5' GCCCGATAATATCACACGAA 3' 








RERG 


F: 5' AGG II II GATCGAGGAAAGTC 3' 


2 


260 


58 




R: 5' CCGACTAACGACGAAAAAAC 3' 








CENPH 


F: 5' ATGGTTGTCGTCGTATTTTC 3' 


3 


320 






R: 5' ACTAAATAAACGCCGAAAAAA 3' 








GSTM5 


F: 5' GTTTTTGGGAATTCGGTAGC 3' 


2 


356 


60 




R: 5' CGAAACCTTACTACGACCCC 3' 








Abbreviations: F, forward primer; R, reverse primer. 
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[3]. Genomic DNA (200 ng) was modified by using the Sigma 
Imprint™ DNA modification kit (Sigma- Aldrich, St. Louis, 
MO, USA) according to the manufacturer's instructions. For 
the 5 genes found identified in the IPA analysis, the DNA 
sequences flanking the transcription start sites (-1000 to +300 
bp), which usually contain a promoter region, were obtained 
from the lUumina Infinium database. Methylation-specific 
PGR primers were designed using the Methyl Primer Express 
Software V.1.0 (Applied Biosystems, Foster Gity, GA, USA). 
The modified DNA from each gene was amplified by 2 differ- 
ent primer pairs specific to unmethylated or methylated 
sequences. The primer sequences used to amplify the GpG 
regions of each gene are listed in Table 1. PGR conditions 
were optimized for accurate and efficient amplification. 
Generally, amplification was carried out in a 25-|liL reaction 
mixture that contained 2 |liL bisulfite -modified DNA, 0.4 
|LiM of each forward and reverse primer, 10 mM Tris-HGl 
(pH 8.3), 50 mM KGl, 1.5 mM MgGl2, 200 |liM dNTPs, and 
2.5 units of Taq DNA polymerase (Qiagen). PGR conditions 
included initial activation and denaturation at 95°G for 15 
min, 35 cycles of amplification (95°G for 45 sec, annealing 
temperature for 45 sec, and 72°G for 1 min), and a final 
extension at 72°G for 10 min. The amplified products were 
resolved on a 2% agarose gel and stained with ethidium 
bromide to enable visualization by ultraviolet light. 



RESULTS 



1 . Patient characteristics 

Genomic DNA was extracted from the bone marrow aspi- 
rate samples for initial diagnosis. Table 2 lists patient charac- 



teristics (N=20) in this study. The median age was 58.5 years 
(range, 16-75), and 13 patients (54.2%) were men. Gharacter- 
istics like sex and age did not differ between the RAEB 
and RGMD groups. 

2. Different methylation profiles between RAEB and RCMD 

According to hierarchical clustering, RGMD patient MDS 
13 had to be excluded from statistical analysis because of 
the extreme value (Fig. 1). In addition, hierarchical analysis 
revealed no meaningful hypermethylated gene cluster pat- 
terns (Fig. 2). The Mann- Whitney and Kmskal-Wallis tests 
were used to identify 1,311 and 1,381 statistically significant 
genes (P<0.05), respectively. A bootstrap t-test and ANOVA 
analysis revealed 1,527 and 821 statistically significant genes 
(P<0.05), respectively. The investigators identified 320 
overlapping genes that displayed significantly different DNA 
methylation patterns between RAEB and RGMD samples. 
Among the 320 genes, RAEB genes exhibited a p-value great- 
er than 0.5, whereas RGMD genes exhibited a p-value less 
than 0.5. Based on any of the 4 statistical models, 69 genes 
exhibited significant differences in the methylation status. 

3. IPA-based selection of candidate genes 

The 69 significant genes were mapped to 5 different net- 
works according to IPA scores (Table 3) as foUows: network 
1: gene expression, cancer, and cell cycle; network 2: organ- 
ismal injury and abnormalities, cellular function and main- 
tenance, and cell death; network 3: cardiovascular disease, 
cancer, and cellular movement; network 4: amino acid me- 
tabolism, post-translational modification, and small molecule 
biochemistry; and network 5: cancer, genetic disorder, and 
hematological disease. Network 1 had a top score of 34 and 



Table 2. Patient characteristics (N=20). 



Sample 


Gender 


Age (yr) 


WHO classification 


BM blast (%) 


IPSS 


WPSS 


MDS6 


Male 


34 


RCMD 


0 


lntermediate-1 


Low 


MDS15 


Female 


35 


RCMD 


0 


lntermediate-1 


Low 


MDS22 


Male 


51 


RCMD 


0 


lntermediate-1 


Low 


MDS7 


Male 


66 


RCMD 


1.0 


lntermediate-1 


Low 


MDS5 


Male 


68 


RCMD 


1.2 


lntermediate-2 


High 


MDS18 


Male 


51 


RCMD 


1.4 


lntermediate-2 


High 


MDS3 


Female 


66 


RCMD 


1.6 


Low 


Low 


MDS8 


Female 


25 


RCMD 


3.6 


lntermediate-1 


Intermediate 


MDS10 


Female 


64 


RCMD 


4.0 


lntermediate-2 


High 


MDS13 


Male 


58 


RCMD 


4.4 


lntermediate-2 


High 


MDS2 


Male 


69 


RCMD 


4.5 


lntermediate-2 


High 


MDS21 


Female 


51 


RAEB-1 


5.8 


lntermediate-2 


Very high 


MDS1 


Male 


62 


RAEB-1 


5.8 


lntermediate-2 


Very high 


MDS17 


Male 


36 


RAEB-1 


6.6 


lntermediate-1 


High 


MDS4 


Male 


47 


RAEB-1 


7.0 


lntermediate-1 


Intermediate 


MDS14 


Male 


67 


RAEB-2 


10.8 


lntermediate-2 


High 


MDS19 


Male 


75 


RAEB-2 


11.0 


lntermediate-2 


High 


MDS12 


Male 


59 


RAEB-2 


11.4 


lntermediate-2 


High 


MDS9 


Female 


16 


RAEB-2 


14.6 


High 


High 


MDS16 


Female 


65 


RAEB-2 


16.6 


lntermediate-2 


High 



Abbreviations: MDS, myelodysplastic syndrome; RCMD, refractory cytopenia with multilineage dysplasia; RAEB, refractory anemia with excess 
blast; IPSS, international prognostic scoring system; BM, bone marrow; WPSS, WHO classification-based prognostic scoring system. 
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Fig. 1 . Cluster pattern for all microarray gene data. Abbreviations: RCMD, refractory cytopenia v^ith multilineage dysplasia; RAEB, refractory anemia 
with excess blast. 




Fig. 2. Hierarchical clustering analysis. No notev\/orthy hypermethylat- 
ed gene cluster patterns v^ere identified. 

is presented in Fig. 3; it includes the associated genes, which 
are denoted in red and white and burrowed to form the 
disease network construct. In total, 16 network 1 genes were 
analyzed (Table 4). These genes included several potential 



links to disease progression in malignancy, and of these, 
5 candidate genes were selected for further evaluation: gluta- 
thione S-transferase mu 5 {GSTM5); BCL2 -interacting killer 
{BIK); centromere protein H {CENPH)\ ras-related, estro- 
gen-regulated growth inhibitor {RERG)\ and angiopoie- 
tin-like protein 2 {ANGPTL2). 

4. Validation of candidate genes by MSP analysis 

To validate the results obtained through methylation mi- 
croarray, we used methylation-enriched PGR to detect meth- 
ylation levels of candidate genes in RCMD and RAEB. The 
5 aforementioned genes were evaluated by MSP analysis. 
Table 1 lists the primer sequences, the optimized annealing 
temperature for each primer pair, and the PGR product 
lengths. The primers were designed to detect the methylated 
and unmethylated DNA sequences of each gene. To ensure 
accurate and efficient amplification of methylated DNA by 
the methylation-specific primers, we optimized the anneal- 
ing temperatures for these primer pairs by using the human 
leukemic-60 cell line as a positive control and peripheral 
blood from healthy volunteers as a negative control. CpG 
lesions of the CENPH gene did not result in methylated 
PGR products following reactions with the specific primer 
pair, as compared to the other 4 genes. 

The MSP methylation results of each gene are shown 
in Fig. 4. RGMD patients were designated numbers 1-10 
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Table 3. Top networks based on I PA analysis. 
Network Associated network functions Score 



1 Gene expression, cancer, cell cycle 34 

2 Organismal injury and abnormalities, cellular function and maintenance, cell death 31 

3 Cardiovascular disease, cancer, cellular movement 29 

4 Amino acid metabolism, post-translational modification, small molecule biochemistry 12 

5 Cancer, genetic disorder, hematological disease 2 




and RAEB patients, 11-19. In the 6^5771/5 promoter, 2 of 
the 10 RCMD samples (20%) and 5 of the 9 RAEB samples 
(55.5%) were shown to harbor different levels of methylation 
sites. In the ANGPTL2 promoter, 1 of the 10 RCMD samples 
(10%) and 4 of the 9 RAEB samples (44.4%) were shown 
to harbor methylation sites. The BIK promoter exhibited 
methylation in 2 of the 9 RAEB samples (22%). The RERG 
promoter exhibited methylation in 3 of the 10 RCMD samples 
(30%) and 2 of the 9 RAEB samples (22%). 



DISCUSSION 



DNA methylation is associated with the transcriptional 
silencing of genes in cell differentiation, apoptosis, growth 
regulation, cell signaling, detoxification, and tumor invasive- 
ness and thus confers a selective growth advantage to malig- 



nant cells [4]. AML and MDS are diseases associated with 
aberrant methylation in a number of genes [5, 6]. Hypermeth- 
ylation of genes involved in cell-cycle control and apoptosis 
is observed commonly in these diseases, particularly 
high-risk MDS. Transcription of the pl5 gene often is si- 
lenced due to abnormal methylation of its promoter region, 
and several studies have indicated that approximately 50% 
of MDS patients exhibit this alteration [7, 8]. Other genes 
frequently affected by hypermethylation in MDS are cancer 
1 {HICI), E-cadherin {CDHI), and estrogen receptor {ER) 
[9]. 

The hypermethylation of promoters from a set of can- 
cer-related genes has been linked to the progression of MDS 
[10]. By examining 1,505 CpG sites in 807 genes using micro- 
array technique and a high-density single nucleotide poly- 
morphism array, Jiang et al. reported that aberrant hyper- 
methylation is associated with disease stage in MDS. The 



Korean J Hematol 2012;47:186-93. 



Different methylation profiles of MDS 



191 



ubiquity, extent, and progression of aberrant hypermethy- 
lation in the disease stage suggests that hypermethylation 
is the primary source of tumor suppressor gene silencing 
and clonal variation in the neoplastic evolution of myeloid 
disease. 

The present study focuses on bone marrow blast in primary 
MDS because increased numbers of bone marrow blasts occur 
during disease progression and are linked to an increased 
risk for transformation to AML. The hypothesis that drove 



Table 4. Statistically different gene methylation in network 1 . 



Gene symbol 


Average p-value 
in RAEB 


Average p-value 
in RCMD 


GSTM5 


0.66136 


0.41475 


BIK 


0.63327 


0.49143 


TNS1 


0.62063 


0.49668 


ABCG4 


0.61621 


0.47086 


CDH15 


0.60316 


0.42345 


FNDC5 


0.57755 


0.4722 6 


AQP1 


0.57256 


0.48058 


CENPH 


0.55901 


0.42925 


RERG 


0.55842 


0.41848 


ANGPTL2 


0.54428 


0.39866 


MEGF10 


0.53940 


0.46604 


EPS8L1 


0.53417 


0.40732 


FGF20 


0.53212 


0.29549 


KIAA1267 


0.51957 


0.47072 


DYDC1 


0.51597 


0.31463 


MALL 


0.50756 


0.24656 


Abbreviations: RCMD, refractory cytopenia with multilineage dys- 
plasia; RAEB, refractory anemia with excess blast. 



the present study was that the mechanism of increased bone 
marrow blast will provide clues to the neoplastic evolution 
of MDS and may be associated with aberrant DNA methyl- 
ation. Based on the WHO classification, 9 patients of the 
RAEB subtype and 11 patients of the RCMD subtype were 
selected and used to determine the different patterns of gene 
methylation status between the 2 subtypes. 

The use of complimentary DNA microarray technology 
in MDS has demonstrated a high level of gene expression 
heterogeneity among MDS patients [11]. However, the 
profiling of DNA methylation in MDS by the microarray 
technique has not been studied yet. The advent of microarray 
technology has provided new opportunities for high-throu- 
ghput study on DNA methylation [12]. The whole ge- 
nome-wide microarray approach has been used to study the 
mechanisms of cancer disease progression and prognosis. 
Thus, the present study investigated DNA methylation in 
MDS by using a commercial whole-genome-wide DNA mi- 
croarray technique. The Infinium Methylation Assay is a 
new technology that differs from the Goldengate Methyl- 
ation Assay in terms of whole-genome-wide content, but 
not cancer/custom content. It still is not practical to use 
genome-wide sequencing of individual tumors to assess the 
entire mutation spectrum. Similarly, research that uses the 
genome-wide DNA methylation assay is in its early stages 
and expensive for use. These considerations underlie the 
decision to limit the assay to 20 MDS patients in this initial 
study. 

Among the 14,495 genes, 69 were hypermethylated in 
RAEB and hypomethylated in RCMD. The DNA methylation 
profile of the 20 samples exhibited a high level of hetero- 
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M 
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11 
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M 
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Fig. 4. Methylation-specific PCR 
results on GSTM5, ANGPTL2, BIK, 
and RERG genes from bone 
marrow samples. Patients were 
designated numbers 1-10 (RCMD) 
and 11-19 (RAEB). The PCR pro- 
ducts were analyzed by electro- 
phoresis on a 2% agarose gel. 
Abbreviations: N, normal peri- 
pheral blood; H, Human leukemic 
HL-60 cells; M, amplified products 
used as primers for the methylated 
sequence; U, amplified products 
used as primers for the unme- 
thylated sequence. 
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geneity. The small sample size limited identification of differ- 
ent patterns between subtypes based on hierarchical cluster- 
ing. Furthermore, this study was conducted on mononuclear 
nonadherent bone marrow cells from routine aspirates. These 
samples represented a non-homogeneous pool of different 
maturation stages, which resulted in heterogeneity. 

It is a daunting task to use microarray studies to discover 
biological links between these genes and disease progression 
and to integrate them within a connected network. The 
present study utilized IPA software to analyze the candidate 
genes further. The IPA-computed score of each network 
indicates the likelihood that input genes in the network 
will be found together due to random chance. The 69 candi- 
date genes analyzed using IPA mapped in 5 networks (see 
Table 3). The top scoring network (network 1) was selected 
for further study. 

Although the functions of the 5 candidate network genes 
{GSTM5, BIK, CENPH, RERG, and ANGPTL2) were not 
determined, previous studies have provided sufficient evi- 
dence that links these genes with malignancy. The methyl- 
ation status of each gene was validated by MPS analysis 
on both subtypes. 

The GSTM5 gene encodes a GST that belongs to the mu 
class of enzymes, which conjugates with glutathione to detox- 
ify electrophilic compounds, including carcinogens, ther- 
apeutic drugs, environmental toxins, and products of oxida- 
tive stress. Genetic variation in mu genes can change an 
individual's susceptibility to carcinogens and toxins and also 
affect the toxicity and efficacy of certain drugs [13]. Express- 
ion of mu class GST has been associated with a poor prognosis 
in acute lymphoblastic leukemia that results from cytotoxic 
drug resistance [14, 15]. Recently, Peng et al. [16] reported 
that DNA hypermethylation regulates the expression of 
mu-class GST in Barrett's adenocarcinoma. In the present 
study, MSP analysis results indicated that the GSTM5 gene 
was methylated more frequently in RAEB patients than in 
RCMD patients (55.5% vs. 20%, respectively). However, in 
comparison with the normal control, slightly methylated 
products were also noted by MSP. Although the possible 
role for GSTM5 in the development of cancer remains un- 
certain, the hypermethylation of GSTM5 may be associated 
with the function of other members of the mu-class GSTM 
family and may change individual susceptibility to carci- 
nogens and toxins, thereby resulting in secondary MDS muta- 
tions that underlie disease progression. 

The ANGPTL2 gene encodes angiopoietin-like protein 2, 
a secreted protein belonging to the angiopoietin family. 
Kikuchi et al. [17] have reported that y47V6!P7Z2 frequently 
is silenced in ovarian cancer cell lines via DNA methylation 
of sequences around the CpG-island that exhibit promoter 
activity. Zhang et al [18] reported that ANGPTL2 and 
ANGPTL3 proteins stimulate ex vivo expansion of hema- 
topoietic stem cells. Here, the ANGPTL2 gene displays a 
slightly higher methylation status in RAEB patients than 
in RCMD patients (44.4% versus 10%, respectively), but 
both subtypes exhibited less DNA methylation. Although 
the physiologic and pathologic functions of ANGPTL2 re- 



main to be discovered, epigenetic silencing by hyper- 
methylation of the ANGPTL2 promoter may lead to a loss 
of ANGPTL2 function and affect hematopoietic stem cells 
in the bone marrow. ANGPTL2 may function differently 
in leukemic and myelodysplastic cells. 

Ectopic expression of BIK induces apoptosis in numerous 
tumor types, including colon, prostate, breast, melanoma, 
and glioma cells [19-21]. BIK interacts with Bcl-xL and Bcl-2 
and triggers apoptosis in a Bak-dependent manner [20]. BIK 
may be a tumor suppressor gene that is downregulated by 
hypermethylation or histone deacetylation of the promoter 
[22, 23]. In this study, a methylated allele of BIKwds observed 
in 2 samples of RAEB, but not in any RCMD samples. In 
other studies, increased apoptosis has been documented with- 
in CD34^ bone marrow cells in patients with MDS [24]. 
Increased apoptosis appears to be the highest in patients 
with refractory anemia and is associated with early stage 
MDS [25, 26]. Although methylation of the BIK gene was 
detected less frequently in RAEB patients, it is possible that 
BIK expression may be downregulated by promoter hyper- 
methylation. Decreased levels of apoptosis in advanced MDS 
may be associated with BIK expression. 

RERG is a novel gene that was first reported in breast 
cancer following microarray analysis [27]. Gender-depend- 
ent expression of RERG in hepatocellular carcinoma and 
its regulation by histone deacetyltransferases (HDACs) have 
been reported previously [28]. In this study, the RERG gene 
was selected as a different candidate gene between RCMD 
and RAEB samples, but no specific pattern was noted upon 
MSP analysis. The Ras family has been studied in MDS and 
leukemogenesis, and some of these mutations may influence 
the course of disease. Therefore, the RERG gene may be 
associated with myeloid disorder and regulated by multiple 
pathways, including the HDAC-related signaling pathway. 

The methylated primers, which were biased to detect me- 
thylated alleles, detected alleles that were 100% hyper- 
methylated. Because MSP analysis allows for only a small 
number of individual CpG dinucleotides to be sampled, a 
negative result does not exclude methylation elsewhere in 
the CpG island [29]. Thus, it was necessary to conduct bisul- 
phate sequencing analysis of the CpG island to determine 
the methylation status of each gene. 

In this study, the whole-genome-wide DNA microarray 
technique and IPA were used to determine the differing 
states of DNA methylation between RAEB and RCMD; 16 
genes were shown to exhibit increased hypermethylation 
in RAEB samples compared to RCMD samples. Among the 
16 genes, GSTM5, ANGPTI2, and ^/iT promoters exhibited 
methylation more frequently in RAEB than in RCMD 
samples. Although the functions of these genes remain un- 
known, their state of DNA methylation may induce epi- 
genetic silencing and contribute to increasing blasts, which 
results in the disease progression of MDS. Additional studies 
focusing on epigenetic silencing by various detection modal- 
ities are required. 
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